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Abstract. We study the weak boundary layer phenomenon of the Navier-Stokes equa- 
tions with generahzed Navier friction boundary conditions, u-n — 0, [S{u)n]^^^^+ Au = 
0, in a bounded domain in R'^ when the viscosity, £ > 0, is small. Here, ^ is a type 
(1, 1) tensor on the boundary: when A = al we obtain Navier boundary conditions, 
and when A is the shape operator we obtain the conditions, u ■ n = (curlw) x n = 0. 
By constructing an explicit corrector, we prove the convergence, as e tends to zero, 
of the Navier-Stokes solutions to the Euler solution. We do this both in the natural 
energy norm with a rate of order e^^* as well as uniformly in time and space with a 
rate of order g^/s-'S near the boundary and e^/'^"^ in the interior, where 5, 5' decrease 
to as the regularity of the initial velocity increases. This work simplifies an earlier 
work of Iftimie and Sueur, [22], as we use a simple and explicit corrector (which is more 
easily implemented in numerical applications). It also improves a result of Masmoudi 
and Rousset, [29], who obtain convergence uniformly in time and space via a method 
that does not yield a convergence rate. 
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1. Introduction 

The flow of an incompressible, const ant- density, constant-viscosity Newtonian fluid is 
described by the Navier-Stokes equations. 



eAw^ + {u^ ■ V)m" + Vp'' = f in ^] X (0, T), 

divu^ = infix(0,T), ^^'^^ 
M^|t=o = Mo in ^• 

The fluid is contained in the three-dimensional bounded domain, fi, with smooth bound- 
ary, r. The parameter, e > is the viscosity and T > is fixed (see Theorem 2.2). The 
equations are to be solved for the velocity of the fluid, m^, and pressure, p^, given the 
forcing function, /, and initial velocity, uq. The regularity of F, /, and uq we assume is 
specified in (1.7), but our emphasis is not on optimal regularity requirements. 
When 5 = 0, we formally obtain the Euler equations, 

— + (mO ■ V)m° + V/ = / in X (0, T), 

divM° = infix(0,T), (1.2) 
M°-n = onrx(0,T), 
u\=Q = Uo in fi, 

where n is the outer unit normal vector on F. 

For the Euler equations, we use the minimal, impermeable boundary conditions, 
■ n = 0, reflecting no entry or exit of fluid from the domain; being a first-order 
equation, these conditions suffice. No-slip boundary conditions, = on F, are those 
most often prescribed for the Navier-Stokes equations. This, of necessity, leads to a 
discrepancy between and vP at the boundary, resulting in boundary layer effects. 
Prandtl [34] was the first to make real progress on analyzing these effects, and much of a 
pragmatic nature has been discovered, but to this day the mathematical understanding 
is woefully inadequate. 

In part because of these difficulties with no-slip boundary conditions, and in part 
because of very real physical applications, researchers have turned to other boundary 
conditions. Of particular interest are boundary conditions variously called Navier fric- 
tion, Navier slip, or simply Navier boundary conditions (other names have been used as 
well). These boundary conditions can be written as 

■ n = 0, [S(u^)n + au\^^^ = on F, (1.3) 

where 

1 /_ .^x /15mi Idui 

2 dXj ) \<i,3<Z 

Here (xi, a;2, X3), (or {x,y,z) in Section 3), denotes the Cartesian coordinates of a point 
X e M^, a is the (positive or negative) friction coefficient, which is independent of e. 
The notation in (1.3) denotes the tangential components of a vector on F. 
In this paper, we use the generalization of (1.3), 

■ n = on F, 



S(..) l{Vu + (V.)') ^ (1^ + ig)^ , fo. . ^ (1.4) 
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of the Navier boundary conditions. Here, ^ is a type (1, 1) tensor on the boundary 
having at least C^-regularity. In coordinates on the boundary, A can be written in 
matrix form as .4 = [aij) Note that lies in the tangent plane, as does Au'^. 

It is easy to see that when A = al, the product of a function a on F and the identity 
tensor, the generalized Navier boundary conditions, (1.5), reduce to the usual Navier 
friction boundary conditions, (1.3). In fact, the analysis using a general A in place of 
al is changed only slightly from using al with a a constant (we say a bit more on this 
in Remark 2.5). 

The primary motivation for generalizing Navier boundary conditions in this manner 
is that when A is the shape operator (Weingarten map) on F, one obtains, as a special 
case, the boundary conditions, 

■ n = (curlw^) x n = 0, (1.6) 

as we show in Appendix B. (This fact is implicit in [3].) Such boundary conditions 
have been studied (in 3D) by several authors, including [2, 3, 47] (and see the references 
therein), [7, 6] for an inhomogeneous version of (1.6), and [4, 5] for related boundary con- 
ditions. In this special case, stronger convergence can be obtained (at least in a channel), 
in large part because vorticity can be controlled near the boundary. Hence, somewhat 
different issues arise, and the bodies of literature studying boundary conditions (1.6) 
and (1.3) are somewhat disjoint. 
We introduce the Hilbert space, 

H = [u^ L'^{nf : div M = and M ■ n = on F} , 

equipped with the usual inner product. Then, letting T > be an arbitrary time 
less than any T appearing in Theorems 2.2 and 2.4, we state our main result: 

Theorem 1.1. Let Fa he the interior tubular neighborhood of VL with width a > 0. 
Assume that 

uoeHn H"'{Sl), f G ^'^.([O, oo); C°°(fi)), F is 6""+^ for m > 5. (1.7) 

Then m^, a solution of the Navier-Stokes equations, (1-1), with Generalized Navier bound- 
ary conditions, (1.5), converges to , the solution of the Euler equations, (1-2), as the 
viscous parameter e tends to zero, in the sense that 

for some T > and for a constant n = k(T, a, mq, /), a = WAWQ^f^Y^y If m > 6 and 
/ = then 

M nir 33 .. 39 

ll" \\L-^{[o,T]xra) - ' II " llL°°([o,T]xQ\r„) - '^^ ' y^-^' 

where now k = k{T, a, m, a, txq, /) ■ 

Because we will only have existence of u"^ when 5 < m < 6 (Theorem 2.1), by 7/^ we 
mean an arbitrary choice of the possibly multiple solutions when we consider the limit 
as e — )■ 0. When m > 6 the solutions are unique as shown by Masmoudi and Rousset 
(see Theorem 2.4), and this arbitrary choice becomes unnecessary. 

Remark 1.2. Standard boundary layer analysis indicates that a linear corrector will 
be of order L°°([0,T] X fl), so an exponent of | rather than | in (1.9) should be 

considered optimal (for C°° initial data). 
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Navier boundary conditions go back to [32], in which Navier first proposed them, and 
to [30], in which Maxwell derived them from the kinetic theory of gases. There has been 
intermittent interest in them since, but revival of active interest in the mathematical 
community working on the vanishing viscosity limit started with the paper of Clopeau, 
Mikelic, and Robert [10], which gives a vanishing viscosity result in two dimensions. 
Also, the work of J-M Coron in [11] on the controllability of the 2D Navier-Stokes 
equations with Navier boundary conditions, which precedes [10], initiated interest in 
these boundary conditions in the PDE control theory community. By now there is a fairly 
substantial mathematical literature on the subject, but the three papers, [21, 22, 29], are 
of particular concern to us here. Both [21] and [29] give existence theorems for solutions 
to (1.1, 1.3), with uniqueness holding for stronger initial data. We quote these results 
in Theorems 2.1 and 2.4. 

Even with Navier boundary conditions there is a discrepancy between and on 
the boundary, so we expect boundary layer effects to occur. As first shown (in 3D) 
by Iftimie and Planas in [21], however, this boundary layer effect is mild enough to 
allow convergence of u"^ to in L°°{0,T] L'^{Q)) without using any artificial function 
correcting the difference u"^ — vP on the boundary. (This result of [21] was for a > 0, 
but the argument is easily modified to allow a to be negative.) Thus, it makes sense to 
refer to the boundary layer as weak. 

Specifically, Iftimie and Planas show in [21] that 

\\u^ - u°||L^(o,T;L2(n)) < Ce^. (1.10) 

Iftimie and Sueur [22] use a corrector, to improve the convergence rate in (1.10) to 
e^/^ in this energy norm. More precisely, they consider an asymptotic expansion of u'^ as 
the sum of vP and where w is a corrector whose main part of its tangential components 
is defined as a solution of a linearized Prandtl-type equation. Using the properties of 
they show that 

u' - (m° + v) is order e in L°°(0, T; L^{^)), order in L^{<d, T; H\n)). 
These bounds with estimates on the corrector v then give 

- AL^iO,T;LHn)) ^ C'^'' - A\L^O,T;mm ^ ^'^^ (LH) 

We, on the other hand, use an asymptotic expansion of u"^ in the form u"^ ~ m° + 
where the main part of the explicitly defined corrector, 6'^, exponentially decays from the 
boundary; see (3.8, 3.9) and (4.24, 4.25). Because our corrector is so simple and explicit 
it can be mored readily used in numerical applications than that of [22], which requires 
the solution of a coupled system of linear equations. Both correctors are linear and 
both can be used to obtain order e^^^ convergence in the vanishing viscosity limit, (1.8), 
but Iftimie and Sueur achieve an order of convergence of e for the corrected difference, 
^£ _~ ^]-^jig Q^j^ corrected difference still gives order e^^'^. The tradeoff is simplicity 
of the corrector versus rate of convergence of the corrected velocity. 

We wish to emphasize that the techniques we employ in this paper differ considerably 
from those of [22]. While the approach in both papers originates in the work of Prandtl 
[34], our approach adheres much more closely to a by now well-established approach to 
boundary layer analysis, which we adapt to treat Navier boundary conditions. In this 
regard our arguments will be more familiar to many researchers, and hence, ultimately, 
we believe, easier to incorporate into the existing understanding of boundary layers 
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as they appear in a variety of physical problems. (For a description of the general 
theory of boundary layer analysis, see, for example, [12, 13, 14, 16, 20, 27, 33, 35, 45]. 
Concerning boundary layer analysis related to the Navier-Stokes equations, we refer 
readers to [15, 17, 18, 19, 22, 24, 38, 39, 42, 43, 44].) 

A key aspect of our corrector is that it is coordinate-independent. This not only gives 
it geometric meaning, it removes the need for a partition of unity to patch together the 
corrector defined in charts throughout the boundary layer. Nonetheless, the corrector 
has a particularly simple form in principal curvature coordinates, which we discuss in 
Section 4.2. (Such coordinates are used in much the same way, though for different 
purposes, in [3].) 

We also, in (1.9), obtain convergence uniformly in time and space of order e^/s-s ^gg^j. 
the boundary and e'^/^~^ in the interior, with 5, 5' decreasing as the regularity of the 
initial velocity is increased, by employing an anisotropic embedding inequality developed 
in Section A. We take great advantage of the regularity result of Masmoudi and Rousset 
[29] (Theorem 2.4) to obtain this convergence. The authors of [29] themselves take a 
similar approach; however, the anisotropic inequality they use requires control on norms 
higher than those in (1.11), and this is only sufficient to obtain boundedness of the 
sequence of solutions to (1.1, 1.3). A compactness argument then gives convergence 
uniformly in time and space, though without a rate of convergence. 

The body of this paper is organized as follows: The existence and uniqueness results 
for solutions to the Navier-Stokes equations and Euler equations that we will need are 
given in Section 2. We give the proof of (1.8), the first part of Theorem 1.1, in Sections 3 
and 4. To avoid the geometrical difficulties of a curved boundary, which obscure the 
key ingredients of the argument, we first prove (1.8) for the case of a three-dimensional 
periodic channel domain. We do this in Section 3. Then, in Section 4, as a generaliza- 
tion of Section 3, we treat the case of a bounded domain in with smooth (curved) 
boundary. In Section 5, we present the (very short) proof of (1.9), the second part of 
Theorem 1.1, which relies on the anisotropic Agmon's inequality, which we establish in 
Theorem A. 2. In Appendix B, we prove that (1.5) reduces to (1.6) when A is the shape 
operator. Finally, Section C contains some standard lemmas which we state without 
proof. 

2. Existence and Uniqueness Theorems 

Thanks to Lemma C.l, by applying the Galerkin method, one can construct solutions 
to (1.1) with (1.5) in the following sense, as shown in [22] (see Remark 2.5): 

Theorem 2.1 (Iftimie, Sueur [22]). Assuming thatuo lies in H and f lies in L^(0, T; L'^(Q)^) 
there exists a weak solution G C^(0,T;if) H L^(0, T; if^(f2)'^) of the Navier Stokes 
equations, (1.1), with the generalized Navier friction boundary conditions, (1-5), in the 
sense that, for any G if n C^([0, T] x H), 

I u'-^dxdt + 2e I [ S{u') ■ S{v) dxdt 
iQ Jn Jo Jn 

+2e I I Au^-vdSdt+ / / {u"^ ■ V)u^ ■ v dx dt = / UQ-v\t=odx. 
Jo Jr Jo Jn Jn 

We have the following well-posedness result for solutions to the Euler equations: 
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Theorem 2.2. Suppose that vP lies in H D H^iVt) fl C^''^{VL), jj, in (0, 1], m > 3 is an 
integer, f lies in C^^{[0, oo); C°°{Q)), and T is of class C™+^. Then for some time T > 
there exists a unique solution, , to (1-2) lying in C^([0, T] x fi)nC([0, T]; H"^{Q)). The 
corresponding pressure, p^, lies in L°°(0, T; and is unique up to an additive 

function of time. 

Proof. Combining Theorem 1 and Theorem 2 part 3 of [26] gives the existence, unique- 
ness, and regularity of u when / = and the boundary is smooth. The proof is straight- 
forward to adapt to smooth forcing, and the strongest restriction on the smoothness of 
the boundary comes through the use of the Leray projector (Lemma 2 of [26]), where, 
however, C^^^-regularity suffices. The regularity of the pressure (as well as the well- 
posedness in Sobolev spaces) is proved in [41, 37]. □ 

When (1.7) holds, by virtue of Theorem 2.2 and Sobolev embedding, for some T > 0, 

M° e Cl{[0,T] X 1]) nC([0,T];/7™(fi)), / e L'^{0,T;H'^-^\n)) for m > 5. (2.1) 

The regularity in (2.1) of the solution is what we require; we do not claim that the 
assumptions in (1.7) are the minimal ones that guarantee such regularity, however. 

In [29], Masmoudi and Rousset obtain the well-posedness result that we state in 
Theorem 2.4 for solutions to (1.1, 1.5) in the conormal Sobolev spaces of Definition 2.3 
(see Remark 2.5). 

Definition 2.3. Let Q he a. rf-dimensional manifold, d > 1, with C'^-boundary, k > 1. 
Viewing vector fields as derivations of C°°{Q), we say that a vector field, X, is tangent 
to dfl if Xf = on dQ whenever / is constant on Q. Let {Zj)jLi be a set of generators 
of vector fields tangent to dQ. (Locally, only d such vector fields are needed, but for 
a global basis, will be greater than d.) For a multiindex, /3, let = Zf^ ■ ■ ■ Z^^ . 
Define 

H^^{Q) = {f e L^Q): Z^f G L\Q) for all |/3| < m} 

with 

\l3\<m 

We say that / is in the space, W^'°°, if 

\(3\<m 

and we define the space E'^ by 

:= {u e HZ{Q)\ Vu e HZ^'m 

with the obvious norm. 

Theorem 2.4 (Masmoudi, Rousset [29]). Letm he an integer satisfying m > 6 andVt he 
a (7™+^ domain. Consider uq G E"^nH such that Vuq G W^f°. Then there exists T > 
such that for all sufficiently small e there exists a unique solution, G C([0, T], i?™), to 
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(1.1, 1.5) with / = and such that ||Vm^||i,oo bounded on [0,T]. Moreover, there exists 



C = C{a) > 0, where a 
sup(||u^ 

[0,T] 

+ e 



; such th(xt 
Hsiin) + l|VM"(t)||j^m-i(f^) + ||Vu"(t)||^^i.^) 

\\VMs)\\l^^-^^n)ds<C. 



(2.2) 



Remark 2.5. Theorems 2.1, 2.2, and 2.4 were proved for a bounded domain, but each 
of the proofs extends easily to a 3D channeL Theorems 2.1 and 2.4 were also proved 
assuming that A = al, where a is a constant, but they easily extend to a general A 
by using a = ||^||(7m(r) in place of a in certain boundary terms, much as we do in 
Sections 3.2 and 4.3. 

3. Channel domain 

In this section, we prove (1.8) for a periodic channel domain in R^. We set Q^o '■= 
X (0, h), and consider solutions to (1.1, 1.5) in a channel domain Q^o- That is. 



dt 



eAu^ + {u^ ■ V)m" + Vp' = /, in x (0, T), 



^ div = 0, in n^o x (0,T), 

and are periodic in x and y directions with periods Li and L2, 



(3.1) 



Here, / and uq, satisfying (1.7), are assumed to be periodic in x and y directions with 
periods Li and L2, respectively. 

Since n = (0, 0, —1) at 2; = and n = (0, 0, 1) ai z = h, we can write the Generalized 
Navier boundary condition, appearing in (1.5) with (1.4), in the form 



n = 

dz 

duf 
dz 



0, at z = 0, h, 



2 
2 



0, i 



1, 2, at z 



0, 



(3.2) 



0, i = 1,2, at 2; = h. 



The corresponding limit problem, (1.2), can be written as 
+ (uO ■ V)nO + VpO = /, in x (0, T), 



dt 



div M° = 0, in x (0,T), 
^ M° and p° are periodic in x and y directions with periods Li and L2, 



(3.3) 



u 







0, at z = 0, /i, 



13 — c»,^ ^ — 

I, u°\t=o = Uo, in r^oo- 
For the sake of convenience, we set 



n := (0,Li) X (0,1.2) X (0,/i), 
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and assume (to simplify the expressions in (3.6)) that 

e < {h/sy. 

To study the boundary layer associated with the Navier-Stokes problem (3.1) with 
the Navier friction boundary conditions (3.2), we expand m*^ as 

~ u° + (3.4) 

where m° is the solution of (3.3) and 6^ is a divergence-free corrector, which will be 
determined below. The main role of is to correct the tangential error related to the 
normal derivative of — on the boundary; see (3.5) below. 



3.1. The corrector. To define a corrector, 9^ = {Of, 9^, 6^1), using the ansatz 6*1 ~ e^^^9^, 
i = 1,2, with respect to the order of e in any Sobolev space, we first devote ourselves to 
find a suitable boundary condition for 9f, i = 1,2. By inserting the expansion (3.4) into 
(3.2)2,3, we find that, for i = 1,2, 



dvl 

dz 

dvl 
dz 



09 



- 2 ^ aijU^j + ^ ^ 2 ^ aij9'j ~ 0, at z = 0, 



dz 

89: 



+ 2 ^ aijU^j + ^ + 2 5Z "^J^i - ^' at ;z = /i. 



For smooth Uij, 1 < i,j < 2 on T, independent of e, we expect that 391 /dz ^ 
2 'Yl,j=i Cdj^j^ i = 1,2. Hence, we use the Neumann boundary condition for 9f 



9^ 



091 
dz 

m 

dz 



Ui,L 



Ui,R 



^ - 2 aijU^jj, at z = for i = 1, 2, 



2 



+ 2 '^ij'^'jj 5 a^t z = h for i = 1,2. 



(3.5) 



In the theory of boundary layer analysis, it is well known that the Neumann type 
boundary condition, (3.5), is useful when treating any weak boundary layer phenomenon. 
More precisely, to improve the convergence given in (1.10), it is sufficient to construct 
a corrector function that fixes the normal derivative of the difference — m° on the 
boundary, instead of the difference itself. 

Toward this end, we first define cutoff functions, ai, cr^, belonging to C°°([0, h]), by 



1, 0<z<h/^ 



aR{z) := aL{h - z) 



(3.6) 



0, h/A<z< h. 

Then we define the tangential component 91, i = 1,2, of the corrector 6'^ = {91, ^2, ^3) as 



9',:=9U. + 9: 



i,Ri 



1,2, 



(3.7) 
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where 



i,L 



d 



= -euix{x,y]t) — ^aL[z)[l-e 

h — z / h — z 

(^i,R = v^Ui^R{x,y;t)crR{z)e ^ - eui^R{x,y;t)aji{z)\l - e ^ 
= -eUi^R{x,y;t)-^^^aR{z)(^l - e"^j|. 
To make 6"^ divergence-free, we must define the normal component 6^ of 6'^ as 

ne ne , ne 

where 



e|^ + ^)(x,,;tK(.)(l-e-7.), 



dx 



dx 
dui^R 



+ 



dy 

dU2,R 

dy 



[x,y;t)aR{z)(l 



e 



(3i 



(3.9) 



(This form of the corrector is as in [24], adapted to Navier boundary conditions.) 

Thanks to (3.6, 3.8), by differentiating (3.7) with respect to the normal variable z, 
one can easily verify that the tangential components, 61,02, satisfy the desired boundary 
condition (3.5). Moreover, from (3.9), we infer that 



9^ = 0, at z = 0, h. 



(3.10) 



3.2. Bounds on the corrector. We introduce the following convenient notation: 



gk 



any differential operator of order k 

with respect to the tangential variables x and y 



k>0. 



We also use the convention that = kt{T,uq, f) is a constant that depends on T, uq, 
and /, but is independent of e and A, and may vary from occurrence to occurrence. 
Letting 



a 



(3.5) through (3.9) give 

< kt{1 + a)e 



Or 



L°°([0,T]xn) 



C'"(r) 



del 



,m>6, (3.11) 



dz 



< HT{l + a), i = 1,2, (3.12) 



L°°([0,T]xn) 



dr 



< kt{1 + a)e. 



del 



dz 



< Kr(l + a)e 



(3.13) 



L°°([0,T]xn) 

We have the following bounds on the corrector 



L°°([0,T]xn) 
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Lemma 3.1. Assume (1-7) holds and that k,l,n>0 are integers either I = 1, k = or 
I = 0, < k < 2. Then the corrector, 9'^, defined by (3.7) through (3.9), satisfies 



Ql+k+nQ: 



dt^dr^dz^ 



< C(l + i 



1,2, 



Ql+k0e 



dt^dr^ 



L°°(0,T;L2(n)) 



L°°(0,T;L2(Q)) 



Ql+k+n+lQe^ 



dt^dr^dz'^^ 



3 n 
2 



L°°(0,T;L2(Q)) 



for a constant, C = C (T, I, k,n,Uo, f) . 



(3.14) 



Proof. The assumptions (1.7) give the regularity of m° in (2.1), and since m > 5, this 
allows k to be at least as large as 2. To prove the lemma, using (3.6) through (3.9), we 
first notice that it is sufficient to verify (3.14) with 9f replaced by 9f j^, 1 < i <3. 
For (3.14)i with 9f^, using (3.8)i, we write 



Ql+kQe 



1 d'-'^^UiT 



— £2 . 



{aL{z)-e^a'L{z))e ^ -e 



Ql+k 



-0-t\Z . 



(3.15) 



Then, by differentiating (3.15) n times in the z variable, and using (3.5, 3.6), we find 

Ql+k+nQ£ 



dt'-dr^dz'^ 



< C(l + a)£:5-2e + C{1 + a)e + e.s.t., 



(3.16) 



where e.s.t. is a function (or a constant) whose norm in all Sobolev spaces H'^ (and thus 
spaces C^) is exponentially small with a bound of the form ci exp(— 02/^^), Ci, C2, 7 > 0, 
for each s. Hence, we find 



gl+k+nQ, 



e 

i,L 



dt^dr'^dz'' 



< C(l + a)£5 



1 _ n 

2 2 



L°°(0,T;L2{n)) 



e ^ dz] +C{l + a)e 



< (setting z' = z/y/e) 



3 n 
■ 2 



e-^'' dz'Y + C(l + a)e 



<C(l + a)e3-t, for/,A;,n>0. 
To prove (3.14)2 with 6^1^, using (3.9)i, we write 



(3.17) 



& 



l+k 



dUi^L ^ dU2,L 



dy 



aL(z)[l-e ^ 



Hence, (3.14)2 with 6^1^ follows by applying exactly the same computations as (3.16, 
3.17), and the proof of Lemma 3.1 is complete. □ 



We define a continuous piecewise linear function, Cl-^)? by 



z, 0< z< h/A, 
h/4:, h/A<z< 3/i/4, 
h- z, 3h/4: <z <h. 



(3.18) 
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Then, using the analog of the proof of Lemma 3.1, one can verify that, i = 1,2, 



Ciz) del 



y/e dz 



Ciz) del 



a/5 dz 



<C(l + a)£2. (3.19) 



L°°{0,T;L2(r2)) 



< C{l + a)e^, 

L°°(0,T;L2(n)) 

3.3. Error analysis. We set the remainder: 

■=u' -u"^ -e'. (3.20) 
Then, using (3.1) through (3.3) with (3.5, 3.10, 3.20), the equations for read 

-eAw' + V{p' ~p^) =£AMO + i?,(r) - J,(u^uO), in Q^o x (0,T), 



dt 

div = 0, in Qoo X (0,r), 

w'^ is periodic in x and y directions with periods Li and L2, 



wl = 


: 0, at z = 0, h. 






dwi 
dz 


2 

- 2 ^ aijw] = 
i=i 


2 

2 ^ aijej, i = 1,2, at z = 0, 


(3.21) 


dwf 
dz 


2 

+ 2 ^ ai.w'j = 


2 

—2 Oij^*^, 2 = 1,2, at z = /i. 






=0 = -e''\t=o, in 







where 



dv 

Re{v) := + ^At>, for any smooth vector field v. 



(3.22) 

J,(M^ M°) := {u' ■ V)u' - (m° ■ V)m°. (3.23) 

We multiply (3.21)i by w'^ , integrate over VL and then, integrate it by parts . As a result, 
after applying the Schwarz and Young inequalities as well, we find: 
d 



^11^ llL2(n) 



L2{n) 



<e'\\/\u''\\^,^^^ + \\Rem\\l,^^^ + 2\\w 



lL2(n) 



+ 2e I [Vw'n) ■ w' dS - 2 / J,{u' ,u°) ■ dx. 

J{z=0,h} Jq 

Thanks to Lemma 3.1 and (3.22) with v replaced by we find that 



ii?.(mii2(f,)<Mi + «V- 



(3.24) 



(3.25) 



On the other hand, by remembering that n = (0, 0, —1) at 2; = and n = (0, 0, 1) at 
z = h , we notice that 



tan 



dwl dw2 
dz ' dz 

dwl '5^2 



dz ^ dz J' 
(using (3.21)5,6) 



at z = 0, 



at z = h. 



(3.26) 



-2(^^aij{w'^ + e]), ^a2j{w'^ + e])^, at z = 0,h. 
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Then, using (3.26), we find that 
2e 



J{z=0,h} 



[W J tan 1 1 2,2 (D 



< 



KTae ||w^||^2(r) + i^t^s ||[6'^]tan| 



L2(r) ll""^ llL2(r) 



< (using Lemma C.2, (3.8), and the Poincare inequahty) 

< e II Vw"||^2(^) + Kra^e Ww'Wl^f^^^ + kt{1 + a^)e^. 
By applying (3.25, 3.27) to (3.24), we obtain 

^ Ik1li2(c) +^l|Vw"llL2(n) 

<KT{l + a^)e^^ + KT{l + a^) ||?i;^||^2(f^) - 2 / J^{u\u^) ■ dx. 



(3.27) 



(3.28) 



To estimate the last term on the right-hand side of (3.28), using (3.20, 3.23), we first 
notice that 

J,(m", m°) = {u' ■ V)w' + {w' ■ V)(m" - w') + (m° ■ V)r + (r ■ V)m° + (^ ■ V)^ . (3.29) 



Then, we write: 



5 

/ Je{u',U^)-w'dx := y^Je, 



(3.30) 



where 



j;i = hu' ■ V)w' ■w'dx = 0, = hw' ■ V){u' - w') ■ dx 



< = /" (r ■ v)u° ■ dx, = [ {u°- v)r ■ dx, 

Jn Jn 

Jl = I {9' ■ V)r ■ dx. 
Jq 

To bound JT"/, using (3.20), we first write 

= I [w' ■ V)m° ■w'dx+ [ (w' ■ V)r ■ w' dx. 
Jn Jn 



(3.31) 



Then 



{W^ ■ V)«" -W^dx < KT ||V«"||^^(^) MUn) < WWU^) 

and, thanks to (3.12, 3.13), 

/ {w' ■V)e' -w' dx < ktWVO' 
Jn 



'V^in) lk'llL2(n) < /«t(1 + a) \\w'\\l2^n) ■ 



Thus, 



\Je\\ < HT{l + a)\\w' 



lL2(n) 



(3.32) 
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Using Lemma 3.1, we bound by 



L°°(f7) 



IL2(Q) 



\W'-\\r'2,^^ < Kt ||6'^llr2m^ WW" 



\L2{n) 



lL2(Q) 
,£||2 



Il2(Q) 



Since -u^ vanishes at z = or h, using the regularity of we bound JT^"^ by 



(3.33) 



3 

lL°°(n) 

3 



d9' 

J_ _|_ J_ 

9x 9?/ 



L2(Q) 



i llL2(n) 



(3.34) 



"3 




Ciz) de^^ 


C(^) 




\fe dz 



\w 



'J llL2(n) 



< (using (3.18, 3.19) and Lemma 3.1) 
Thanks to (3.12, 3.13) and Lemma 3.1, was can bound by 



86' 



dy J 



Wj dx\ 



E 



09' 
dz ^ 



<E{ii^^ 



1IIL°°(Q) 
3 



89'^ 



dx 



\w 



L2(Q) 



illL2(n) 



+ 11^: 



2llL°=(Q) 



9?/ 



\W 



L2(n) 



J llL2(n) 



(3.35) 



^2; 



\w 



L2(n) 



i llL2(n) 



+ XI ll^3llL-(f7) 

< /tT(l + a^)e^ lk1lL2(f^) 
Then, using (3.32) through (3.35), (3.30) gives 



"II < fi:T(l + a^)e^ + ||w"ii' 



L2(Q) • 



/ Js{u',u°) ■ w' dx < kt(1 + 0^)52 + kt(1 + «) ||w 
in 



£||2 

L2(n) 



(3.36) 



Applying (3.36) to (3.28), we obtain 

^ Ik1li2(n) + £ II Vw"||i2(o) < «;t(1 + a*)£^ + f^Al + a ; 11 11 l^u) ■ 
Moreover, using (3.7) through (3.9) and (3.21)7, we observe that 

1 z_ 3 

lk^|t=o|U2(f7) = 1 1 6''' 1 1 1 L2(!^) < kt{1 + a)e^\\e v^||L2(f^) + l.o.t. < kt(1 + «)£*. 
Thanks to the Gronwall inequality, we finally have the bounds on the remainder, w', 

3 1 

II^1Il-(o,T;l2(q)) < n{T,a,uoJ)e^, \\^'\\L^o,T;mm < K{T,a,uoJ)e^. (3.37) 
3.4. Proof of convergence. Using (3.20), we first notice that 

\u' -u°\<\w'\ + \9'\ pointwise in (]oo X (0,T). (3.38) 
Then, using (3.37) and Lemma 3.1, (1.8) follows from (3.38). □ 
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4. Bounded domain 

In this section we consider the Navier-Stokes equations, (1.1, 1.5), and the Euler 
equations, (1.2), in a bounded domain Q in with boundary F having regularity as 
in (1.7). To handle the geometric difficulties of a curved boundary, we must treat Q as 
a manifold with boundary, ffist constructing charts on F = dQ in a special way, as we 
describe below. 

We consider the boundary, F as a submanifold of M^. Then, since F is a compact and 
smooth surface in M^, we construct a system of finitely many charts where each chart is 
a C™-map, ip, from a domain, U, in to a domain, V, in F. More precisely, we choose 
an orthogonal curvilinear system (^') = (.^1,^2) in U so that, for any point xonV(ZT, 
we write 

x = ^{0, r = (ei,6)e^. (4.1) 

Differentiating (4.1) with respect to ^j, z = 1,2, variables, we obtain the covariant basis 
on U and the metric tensor: 

MC)--=Q^^, ^ = 1,2, (4.2) 

and 

{9ijiO)i<ij<2 ■= [di ■ ^j)i<i,,<2 = diag(^i ■ ^1, 92 ■ 92)- (4-3) 
Moreover, we see that the determinant of the metric tensor is strictly positive; 

g{C,') := det(gij) > 0, for all ^' in the closure of U. (4.4) 

For any smooth 2d compact manifold F in M'^, one can construct a system of finitely 
many charts, which satisfy (4.3, 4.4). Hence, the class of domains under consideration in 
this article covers all smooth bounded domains in with boundary having regularity as 
in (1.7). Moreover, as we will see below in Section 4.2, the construction of the corrector 
is independent of our choice of charts. Hence, it is sufficient to restrict our attention to 
a single chart only, since any estimates we develop will apply equally to all of Q. 

We define Fc to be the interior tubular neighborhood of Q with width c for any 
sufficiently small c > 0, and let a > be small enough that T^a is such a tubular 
neighborhood. We can globally define the coordinate ^3 on T^a to be distance from the 
boundary, with positive distances directed inward. 

We fix the orientation of ^' variables on V so that 

n(0 := -P^MiO), (4.5) 

where ti{C,') is the unit outer normal vector on V. Then, letting U = Ux[0, 3a), we define 
a chart tp: U (giving what are sometimes called boundary normal coordinates): 

X = = ^(eO - e3n(e'), X = (xi, X2, X3) E Fsa- (4.6) 

By differentiating ip in ^ variables and using (4.2), we define the covariant basis of the 
curvilinear system ^: 

9M)=Me)-^3^p,^ = l,'2, 9s{^) = -n{a (4.7) 
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hence, from (4.5, 4.7), we see that the covariant basis satisfies the right-hand rule. 

One important observation here is that the orthogonahty, on V, of g^, i = 1,2, does 
not imply the orthogonality, in V, of g^, i = 1,2, 3. To see this, we first notice that 

—— = (linear combination of g^ and gig)) i = 1,2. (4.8) 

Thus, g^ ■ g^ = foi i = 1,2, but g^ ■ g2 ^ in general. Consequently, the metric tensor 
i9iji$,))i<i,j<3 ■■= (9i ■9j)i<ij<3 satisfies: 

On = i = 3 and j = 1,2, or z = 1, 2 and ?' = 3, 
' (4.9) 

933 = 1- 

Moreover, thanks to (4.4), by choosing the thickness 3a > of the tubular neighborhood 
T^a small enough, we see that 

g{^) := det{gij)i<ij<3 > for all ^ in the closure of f/ = f/ x (0,3a); (4.10) 

The function, := g^^'^, is the magnitude of the Jacobian determinant of the chart, ijj. 

The matrix of the contravariant metric components are defined in the closure of U as 
well: 

I ( 922 -9i2 \ 
i9'')l<^,J<3 = i9ij)l<i,j<3 = - -912 9u . (4.11) 

^ \ 1 J 

We introduce the normalized covariant vectors: 

ei = 1 < 2 < 3. (4.12) 
\9i\ ^ ' 

Then, for a vector valued function F, defined on U , in the form 

3 
i=l 

one can classically express the divergence operator acting on F in the ^ variable (see [9] 
or [25]) as 



V9~^9^i\V9ii J \f9 9^3 
We write the Laplacian of F as 

AF = J2{s'F + CF, + ^)e,, (4.14) 

where 



_ , linear combination of tangential derivatives 
of Fj, 1 < j < 3, in up to order 2 



DFi = ^proportional to 



m 

0^3 



(4.15) 
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Remark 4.1. Note that the coefficients of and D , 1 < i < 3 in (4-15), are multiples 
of y/g , \/9ii ' ^I\f9iii = 1) 2, gi2, g2i, and their derivatives. Thanks to (4-10), 

all these quantities are well-defined because of the regularity assumed for T in (1.7). 

Remark 4.2. Thanks to (4-9), we notice that the tangential directions are perpendicular 
to the normal direction in the tubular neighborhood T^a- Indeed this property enables us 
to obtain the expression of Laplacian as (4-H, 4-15), which is essentially the same as 
for the case of orthogonal curvilinear system. The explicit expression of Laplacian in 
orthogonal system appears in, e.g., [15]. 

For smooth vector fields F, G : t/ — )■ M^, we consider VpG, the covariant derivative 
of G in the direction F, which gives F ■ VG in the Cartesian coordinate system. More 
precisely, we consider the smooth functions F and G in the form 

3 3 

F = 5^F,e„ G = J2G^ei. 

1=1 i=l 

Then, one can write VpG in the $, variable. 



^ ^ {v.{f„F. ^ §) + + Q.(F : G)}e., (4.16) 



where 

dG dG \ ( P'^^'^^'^^ °f ^ linear combination of \ 
Vi[Fi,F2 : = the tangential components Fi, F2 of F, , 

^■^ \ and sum of tangential derivatives of Gi J k^-^') 

Qi{F : G) = (linear combination of the products FjGk, 1 < j, ^ < 3 ). 

Remark 4.3. Qi{F : G), 1 <i <3, are related to the Christoffel symbols of the second 
kind, which comes from the twisting effects of the curvilinear system ^. For the case of 
an orthogonal system, the explicit expression of (4. 16) is given in Appendix 2 of [1]. 

Using the expression of contravariant components of the strain rate tensor, and by 
remembering, from (4.3, 4.7, 4.9), that the covariant basis (and hence the normalized 
covariant basis) is triply orthogonal on F, we write the generalized Navier boundary 
conditions, (1.5), for F = J2i=i ^i^i 
F3 = 0, at ^3 = 0, 

1 dF / dF \ ^ 

~oTTT^ + -^4^*' + y^"*i^i = 0' direction at ^3 = 0, i = 1,2, 

(4.18) 

where 

jy^ ( p dF-i\ _ / linear combination of the tangential component Fi and \ ,^ 
V d^i / V the derivative in of the normal component -F3 j- \ ■ ) 

Remark 4.4. Thanks to (1.7, 4-10), the coefficients of M.i{F), i = 1,2, are well-defined. 
Concerning an orthogonal system, the explicit expression of Aii{F), i = 1,2, appears on 
p. 115 of [31]. 



BOUNDARY LAYER WITH GENERALIZED NAVIER BOUNDARY CONDITIONS 



17 



4.1. The corrector. In defining the corrector, 9'^ , we parallel the strategy we used in 
Section 3 for a periodic channel domain as closely as possible, employing an asymptotic 
expansion, ~ m° + 6*^, as in (3.4), but adapting the corrector to the curved boundary. 
With the unit vectors, e^, defined as in (4.12), on [/, 0^ can be written 



i=l 



(4.20) 



The tangential components 6'?, i = 1,2, will be constructed to correct the tangential 
discrepancy in the boundary conditions related to the normal derivative of — on 
the boundary. Then the normal component 6^ will be deduced from the divergence-free 
condition on 6"^. 

To define the corrector appearing in (4.20), since -u" ■ n = on F, we first set 

M = 2([S(M°)Ha„ + ^M0), 

defined on all of F, and write, in coordinates. 



u 



(4.21) 



i=l 



Then, we insert the expansion ~ m° + 6^ into the generalized Navier boundary condi- 
tions, (1.5), and, thanks to (4.18)2, we find that, for i = 1,2, 

Using (4.19), we expect that 861/ d^^ > Mi{ei , 861/ d^i) E?=i "u^^ ^ = 1,2, for 
smooth aij, I < i,j < 2, independent of e. Hence, we use the Neumann boundary 



condition for 9f, 



891 

8^3 



6=0 



u,i^';t), z = 1,2. 



(4.22) 



We can now model the corrector after the fiat-space version in (3.8). We define a 
smooth cutoff function, cr{C,3), with 



^(6 



1, 0<6<a, 
0, 6>2a. 



(4.23) 



Letting 



7i 



/gi 



6=0 



we define the tangential components of the corrector by 



OKe,t) ■.= -e- 



1,2. 



(4.24) 
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It follows from (4.13) that 



96 



-4{E|-(-)K':'M«(i-*)}.^^. 

Thus, we can easily ensure that 6^ is divergence-free by letting 

Omt) :=^^{E|-(7.^^.)(r;t)}^(e3) (l-e-^) . (4.25) 

It is easy to see that 9^ vanishes at ^3 = 0, and by differentiating (4.24) and us- 
ing (4.23), we see that each tangential component 6'f, i = 1,2, satisfies the boundary 
condition (4.22) to within order ^/e^. 

u.{i'-t)-^eE{i'-t), (4.26) 



where 



^3=0 



E{i'-^t)=y,{0^{^) ^ u.{i';t). (4.27) 



Due to the presence of a in (4.24, 4.25), we also have 



= 0, 1 < i,i < 3, A; > 0. (4.28) 

6>2a 



Remark 4.5. For the case of a flat boundary F, one can choose curvilinear coordinates 
with the metric tensor (fi'ij)i<jj<3, defined in (4.9), as the identity matrix /sxs, and 
hence y^, and 7^, i = 1,2, appearing in (4.24, 4.25, 4.27), are equal to 1. This 
implies that the expression of the corrector defined by (4.24, 4.25) are identical to (3.8, 
3.9) in a channel domain where the error E in (4.27) is now equal to 0. 

Remark 4.6. The form of our corrector (4.24, 4.25) is similar to the background flow 
in Lemma 1 of [46]. 

4.2. The corrector in principal curvature coordinates. In this section, we express 
the corrector in particularly convenient and geometrically meaningful coordinates called 

principal curvature coordinates. 

We define an umbilical point of F to be a point at which the principal curvatures, 
Ki and K2, are equal (this also includes what some authors refer to as a planar point, 
where both curvatures vanish). By Lemma 3.6.6 of [25], in some neighborhood of any 
non-umbilical point there exists a chart in which the metric tensor of (4.2) is diagonal (as 
is the second fundamental form) and the coordinate lines are parallel to the principal 
directions at each point. Such a chart is also called a principal curvature coordinate 
system. 

For now, we assume that we are working in such a chart, ipp-. Up V '^T , 

^ = ^piv'), v' = ivi, V2) e Up. 
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The corresponding covariant basis and metric tensor are 

= ^> ^ = 1,2, (gy(r/'))i<ij<2 = {Qi ■ qj)i<i,j<2 = diag(qri ■ q^, ■ ^s)- (4-29) 



Using (4.6) with ^' and ip replaced by rj' and ipp, we define a chart ipp from Up = 
Up X (0, 3a) into T^a by 

X = M'n) = Mv') - ^3n{r]'), 7] = (V, ^3) G Up. 

As before, ^3 is the distance from the boundary, which we note does not depend upon 
the choice of the boundary chart. 

In the principal curvature coordinate system on Up, the unit outer normal vector n 
satisfies 

dn , 

TT- = f^iivJQi, « = 1,2. 
or]i 

Hence, differentiating ipp in the r] variables gives the covariant basis of the coordinate 
system r}, 



Qiil) = (1 - ^^^{v')^3)q^{v'), ^ = 1,2, q^^{r)) = -n{r]'). 
Using (4.29, 4.30), the metric tensor {(lij)i<i,j<3 is written in the form 



ij / l<i,j <3 








;i-'^2(V)e3)'g22 
1 



(4.30) 



(4.31) 



with its determinant, q{r]), bounded away from zero. This is guaranteed by simply 
choosing the thickness, 3a > 0, of the tubular neighborhood small enough. 

It is easy to see that the coordinate system, ry, derived from the principal curvature 
coordinate system, satisfies (4.4, 4.9, 4.10). Hence we use the expression of the corrector 
9'^ , (4.20, 4.24, 4.25), in rj coordinates and write 



r = + ^^63, 



where 



^--4(^«^'(-^-*))E{ 



qu 



;l^..(C,)(:-e-*)-L{g|-(^ 

for ei = qj/lgj, i = 1, 2, and 63 = -n. 
Using (4.31), it is easy to see that 



^3=0 



?3=0 



(4.32) 



6=0 



1 



for i = 1,2 and j = 3 — i. 



Then, combining (4.32, 4.33), we find 

= (^fe) (1 - e-^)) Miri)u{r^'; t). 



(4.33) 



(4.34) 
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where M is a smooth type (1, 1) tensor defined, in our coordinates, by 

2 ^ 2 

M(r;)F:=^ ^_^ , . m, F = ^fe. (4.35) 

i=l 3 i\ I A3 -^^ 

On the boundary, the divergence operator is 
Then by (4.31), 

= (1 - /«i(V)6)(l - l^2W)^3)Vt, 
so we can write 6^ in (4.32) as 

1 - e _ 
(1 - Ki(r/')6)(1 - i^2{v )^3) 

Although we assumed in its derivation that we were near a non-umbilical point so 
that we could construct a principal curvature coordinate system, our expression for 6^ 
is perfectly vahd at an umbihcal point, where we simply have Ki = K2 thanks to the 
smoothness of the curvatures in the tangential variables. 

Finally, a straightforward but lengthy calculation, which we omit, shows that (4.24, 
4.25) transforms to (4.34, 4.36) under the change of variables from to ipp, showing 
that our corrector in the form (4.24, 4.25) is coordinate- independent. (This is perhaps 
not immediately obvious, because (4.24, 4.25) involve the metrics both on the boundary 
and in the tubular neighborhood.) 

Remark 4.7. For most smooth, bounded domains in M.^, principal curvature coordinate 
systems can be constructed in the neighborhood of all but at most isolated points; in 
fact, having only isolated umbilical points is (in some sense) generic. And, for instance, 
a sphere, while it consists only of umbilical points, can be covered by two charts, both of 
which use principal curvature coordinates (essentially, spherical coordinates). When such 
coordinates suit the boundary of a domain, the expression for the tangential corrector 
in (4.34, 4.35) is both simpler to calculate and simpler to interpret than the expression 
in (4.24). In such coordinates, the expressions for the differential operators, such as 
div, curl. A, can also be written more simply. Though we used principal curvature 
coordinates in this section to prove that our corrector is independent of the choice of 
charts, we cannot restrict ourselves to such coordinates in the rest of our analysis, as 
that would put constraints on the geometry of the domains that we would be able to 
treat. 

Remark 4.8. As a special example of a smooth bounded domain in M^, consider a solid 
torus, which has no umbilical points on its boundary. Hence, we need only one principal 
curvature coordinate system, if we allow it to be periodic in the tangential variables. In 
this sense, the solid torus is the simplest smooth bounded domain to work with in M^. 
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4.3. Bounds on the corrector. We follow the convention described at the beginning 
of Section 3.2, though now the tangential variables are and and as in (3.11), we 
let a = 



d9f 



dr 



(-,„^P^ , m > 6. Then, from (4.24, 4.25), we infer that 

d9f 



< Kr(l + ce)E 



L°°([0,T]x(7) 



< Kril + a), i = 1,2, (4.37) 



L°°([0,T]x(7) 



and 



del 



dr 



< kt{1 + (y)e, 



L^{[0,T]xU) 



del 



0^3 



< kt(1 + a)e2. 



(4.38) 



L°°([0,T]x(7) 



We now state the estimates on the corrector in the lemma below, which we omit the 
proof as it is essentially the same as that of Lemma 3.1 because of (4.10). 

Lemma 4.9. Assume (1-7) holds and that k,l,n>0 are integers either I = 1, k = or 
I = 0, < k < 2. Then the corrector, 6'^ , defined by (4-24, 4-^5), satisfies 



dt^dr^d^. 



<C\V\-^{l + a)e-i-2, i = l,2, 



Ql + kQS 



dt^dr^ 



L°°{0,T;L'^(U)) 

< C\V\^l + a)e, 



L°°(0,T;L2(;7)) 



dt^dT'^dQ^^ 



< C|l^|5(l + a)£3-t, 



L°°(0,T;L2{(7)) 



for C = C(T,l, k,n,uo, f) > 0, independent of e, A and the measure \ V\ ofV = ip{U). 

In addition to the estimates in Lemma 4.9, as for the case of the channel domain, one 
can easily verify that the corrector 9^ satisfies that, i = 1,2, 



^3 def 



< C\V\^{l + a)e^, 



L°°(0,T;L2(C/)) 



4.4. Error analysis. We set the remainder: 



^3 dO: 



L°°{0,T;L^{U)) 



W 



U 



< C\V\^{l + a)e-^. 

(4.39) 



(4.40) 



Then, using (1.1, 1.2, 1.5, 4.40) and the fact that 6*1 = on V, the equations for read 
( dw^ 



dt 



Je{u',u^), mQx (0,T), 



div = 0, in X (0,T), 
■ n = 0, on r, 

[S(w^)n]^^^ + Aw' = - [S{u^ + 0')n]^^^ - A{u^ + 9'), on T, 
w''\t=o = ~0'\t=o, in Q, 
where -Re(-) and Je{-, ■) are defined by (3.22) and (3.23). 



(4.41) 
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We multiply the equation (4.41)i by w^, integrate it over VL and then use Lemma C.l. 
After applying the Schwarz and Young inequalities to the right-hand side of the resulting 
equation, we find: 

d .. ,,,9 



-r \\w 



dt 



lL2{n) 



+ 4e\\S(w' 



lL2(n) 



+ \\RemrL^n) + '^\\^'\\LHn) 



2 / JJu'u^)-w'dx. 



dS 



(4.42) 



To go further, using the Korn inequality, we first notice that 



„e||2 



,e||2 
llL2(n) 



} > /«s||V 



lL2(!^) - n.S\ II V ay ||^2(f^) "T || uy ||^2(j7) / - n^S || v ti? ||^2(n) 

for a constant, Kg, depending on the domain, but independent of e and a. 

Restricted to the range, V , of any chart, ip, we find, using (3.22) with v replaced by 
e^" and (4.14, 4.15) for each 6^ that 

3 



(4.43) 



\Re{e^)\\i.iy) < 



89' 



dt 



< kt{1 + a^)e2, 

where we also used Remark 4.1 and Lemma 4.9. Since we have a finite number of charts 
on r^a, and since 6'^ is supported in T^a by (4.28), the same estimate holds on Q; namely. 



d^ef 



1=1 



lL2(n) 



< k,t(1 + a )e^. 



To estimate the fourth term in the right-hand side of (4.42), we write 
4e j (^Aw' + [S(m° + r )n] + ^(m° + r )) ■ w' dS 



< KT£:a||w||^2(r) + kt£ 



[S(M° + r)n]^^^ + ^(M° + r 



(4.44) 



(4.45) 



L2(r) 



FllL2(r)- 



On each ^ C F, the rage of the boundary chart, ip, using (4.18, 4.19, 4.21, 4.26), we 
have 

'[S(m° + r)n]^^^ + ^(m° + r)' 



V 



E 

■t=i 

2 

E 

i=l 



U 106' 

-Uj — 

2 



86' 



6=0 



2 I 



where Ui and E are defined by (4.21, 4.27). Since this bound holds for all charts, using 
Remark 4.4 and (4.19, 4.37, 4.38), we find that 



[S{u' + 6')n]^^^ + Aiu' + 6' 



L2(r) 



< kt{1 + c?)e-i. 



(4.46) 
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Thanks to (4.43, 4.44, 4.45, 4.46), applying Lemma C.2 and the Poincare inequahty, 
(4.42) yields that 



II el|2 , r, IIVT £l|2 



< kt{1 + a'^)e2 + Kril + a^) ||w^|Il2(q) 



-2 / J,{u',u^) -w'dx. 
Jn 

To estimate the last term of (4.47), using (3.29), we write 



(4.47) 



Js{u',u'^) -w'dx := ^ 
i=i 



■J 



where JT"/, 1 < j < 5, are given by (3.31). Due to (4.28, 4.37, 4.38) and Lemma 4.9, 
one can easily verifies that J'^ , j = 2, 3, satisfies the same estimate, appearing in (3.32, 
3.33), as for the case of a channel domain. That is. 



(n) ■ 



To bound J'^ and J'^, it is sufficient to work in a single chart, ip: U ^V, as there a 
finite number, N, of them, which just introduces the constant, N. 
For J'^, using (3.3I4, 4.28), we write 



\Je\ < ||(m°- V)r 



\w 



Il2(o) 



< (using (4.16, 4.17) with F, G replaced by O'^, respectively) 



L°°{n) 



i=l 



M° ■ 63 



6 



lL2(C/) 
3 



i^°°(r2a) i = l 



^3 oei 



dei 









IL2(Q) 



LHU) 



< (using (4.39), Lemma 4.9 and the regularity of with {u^ ■ 63)1^3=0 = 0) 
For jTjf, using (4.16, 4.17) with G and F replaced by 6*^, and using (4.28), we find 



\LHV) 
3 



\w 



i=l j=l 



del 






L2{U) ^ 



\w 



lL2(f1) 



3\\l°°{U) 



i=l 



d9' 



d^3 



\w 



lL2(n) 



< (using (4.37, 4.38) and Lemma 4.9) 



< kt(1 + a )£4 \\w^ 



<KT{l + Q')e-^ + \\w%l,^^y 
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Using these bounds on Jl, 1 < i < 5, (4.47) becomes 
Moreover, using (4.24, 4.25, 4.4I5), we see that 

1 ^3 3 

lk^|t=o||L2(f7) = ||6'^|t=o|lL2(r2^) < kt(1 + ||e"v^||i2(r2„) + /.o.t. < /tr(l + a)e4. 
Thanks to the Gronwall inequahty, we finally have the estimates, 

3 1 

lk1lL-(0,T;L2(Q)) < f^iT,a,UoJ)ei, |k1lL2(o,T;Hi(f7)) < ^(T, «, Mq, • (4.48) 

4.5. Proof of convergence. Using (4.40), we first notice that 

lu" -u^\ < \w^ + 1^1, (pointwise in x (0,T)). 
Then, using (4.48) and Lemma 4.9, (1.8) follows. □ 

5. Uniform convergence 

With the estimates we now have, the proof of (1.9) is quite simple. 

Because we assume that m > 6 in (1.7), mq G fl H and (by Sobolev embedding) 
V^o G W^co°°- Hence, both (2.1) and the hypotheses for Theorem 2.4 hold^, so we can 
use (2.1, 2.2) to conclude that 

Then using (1.8), Theorem A. 2, and Remark A. 3, (1.9) follows. 

Remark 5.1. Since also — lies in L°°{[0,T]; W^'°°) by Theorem 2.4, we could use 
the Gagliardo-Nirenberg interpolation inequality (see, for instance, p. 314 of [8]), 

2 ,, „ ,, 3 



5 



L°° — II llL2 11"- ||vyi,°o ' 

to give IIm*^ — M°||j;^oo < Ce^^^'^. This is the same rate that is obtained for m = 6 in (1.9); 
since, however, we require that m > 6, (1-9) always gives a better rate than this. 

Appendix A. An anisotropic Agmon's inequality 

In this section we develop a version of Agmon's inequality in d dimensions, c? = 2 or 3, 
that is suitable for applying to anisotropic problems in which there is more control over 
tangential (horizontal) derivatives than over normal (vertical) derivatives. 

We use the notation A B to mean that A < CB for some constant, C, which may 
depend upon the geometry of an underlying domain but not upon anything else. If C 
depends on some parameter, m, then we write A <^rn B. 

Our starting point is the following simple lemma: 



The hypotheses in Theorem 1.1 are not the minimal ones insuring this. 
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Lemma A.l. Let U be a bounded domain in M*^, d = 1,2,3. For any f in H^{U), 
k> d. 



1-- - 
L°°{U) II/IIl2(C/) ^/'^ = 2, 

L°°((7) 11/ II L2(C/) ll/llHfc(C/) if d='i. 

Proof. Combine the ID Agmons' inequality, WfW^oo^^u) ^ ll/llK^cc/) 11/11 

mons' inequality, \\f\\Lo.(u)<. II/IIl2^(c/) \\f\\%l{uv 3D Agmon's inequality, ||/||ioo(t;) < 

II/IIl^V) II/IIh^V) with the Sobolev interpolation inequality, \\f\\HJ{u) <fe II/IIl2(^) Wf^Hki^uy 
0<j<k. □ 

Theorem A. 2. Let Q be a bounded domain in with C"^^^ -boundary, m > 3, and let 
Ta be the tubular neighborhood of fixed width a > interior to Q. Suppose that f and 
V/ lie in the space H^{Q) of Definition 2.3. Then 



1 - ^ i 



L°°(ra) '^rn,a 11/1112(^7 ll/llH^{n) 
L°°(n\ra) ll/llL2(n) Il/lllf™(f7)- 



ll/llL2(C) + l|V/||Hs(n) 



Proof. We define the chart, tp, as in the beginning of Section 4. In this chart, we can 
define a local conormal basis, (Ai, A2, A3), by Xif{x) = di{f o ip)(^ip~^(^x)), i = 1,2 and 

X^f[x) = ^^,-^i ds{f o ^p){'^jJ~^{x)). We will have need, however, only for Ai and A2. 

It suffices to assume that / lies in C°°{Q). Restricting ourselves to the one chart, ip, 
by Lemma A.l, we have, for any ^ = {C,i,C,2,C,3) in U, 

f o m «m 11/ o ^i-, -, ^3)11^(1) 11/ o -, mkiuo) ■ 

Applying Lemma A.l again, this time in ID with k = 1, gives 
\\foi^i;;mlHUo) = [ f°m.i2.izfdi[di', 

J Uo 

« / \\f om,C2,-)\\L^iO,a)\\f °m.C2,-)\\HH0,a) deide2 



< 



Uo 

I tl M|2 



\ 1 

2 / r „ \ 2 



ii/°^(ei,e^,-)ii22(o,„)C^eweU / ii/°^(e;,e^,-)iiHno,a) ^ewe^ 



^ifomi'd^] ( i_ [\fom\' + \vfm)-d3m\"]d^ 



1 1 

2 / /■ „ „ \ 2 



«a( l^\fomwm)\d^) [Jj\fomr + Nfmw]m)\d^ 



iL2(v) 11^ \\m{v) ■ 

The second ^ followed because the magnitude of the Jacobian determinant, J, is 
bounded away from zero and d^ip is bounded above. Because there are a finite number 



of charts on F^, the bounds are uniform over F^. 
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Similarly, for all multiindices, a = (ai,a2,0) with |a| < m, applying Lemma A.l 
with k = 1 gives 



i^"(/°^)(-,-,e3) 



[\D'^(fo^)f + \dsD'^{fo^l;)f] 



= ll^°/llL2(y) ||X°/||^i(^) . 

But this is true for all I a I < m, so 



V Ju 



< 



1/2 



ll/llL2(n) + l|V/||^™(V) 



1/2 



where we can use the conormal Sobolev space since the only derivative in the normal 
direction occurs in V/ itself. Also, because T is (7™-+^, x(j can be chosen to be C"^~^^(Ta) 
and hence f oip has sufficient smoothness. 
Combining these bounds we have. 



< 



1 i_ 

2 2 m 

LHV) 



1 1_ J_ 

2 2m II f II 2m 



1 i_ 

2 2m 

LHV) 



Summing over all the V gives 



1 

2m 



1 1_ 

2 2m 



L2(n) + ll^/llHs(y) 



1 
2 m 



With W = ^l\Ta, Lemma A.l gives 



1 

2 m 



L2 



2 m 

L2(W^) 



3 
2m 



1 3- 

2 m 

L2(n) 



3 
2m 



Combining these last two inequalities completes the proof. 



□ 



Remark A. 3. It is easy to see that Theorem A. 2 holds as well for a channel domain. 



Remark A. 4. When we apply Theorem A. 2 in Section 5 we have full control on the 
tangential derivatives but can control only one derivative in the normal direction. This 
made the proof of Theorem A. 2 quite simple, as we could apply Lemma A.l in 2D to deal 
with both horizontal derivatives isotropically then use Lemma A.l in ID to deal with 
the single normal derivative. Had we needed to deal with each variable anisotropically 
we would have applied Lemma A.l in ID three successive times. 
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It is worth comparing the inequahty in Theorem A. 2 with the 3D Agmon's anisotropic 
inequahty in Proposition 2.2 of [40], which can be written, for any / in H'^{Q), as 

ll/llL°°(n) ^ 

n(ii/ii^^(n)+ii^./iiL^(f.)+ii^.^./ii L2(n) I 

i=i (A.l) 

« WfwUn) {\\f\ym + llv/IL.(n) + IIA/ll) ■ 

In [40], the authors have some control of the Laplacian but not (directly) of the full 
norm. This inequality would not work for us, however, as it includes d^f. 

Both Theorem A. 2 and the inequality in (A.l) are descendants in spirit of Solonnikov's 
Theorem 4 of [36]. The proof in (A.l) uses, in part, Solonnikov's approach. The approach 
we have taken is, however, more elementary and direct than that of [36]. 

Another type of anisotropic inequality that is not a descendant of Solonnikov's theo- 
rem (and is not of Agmon type) is the anisotropic embedding inequality of [24] (Corollary 
7.3), which originated in Remark 4.2 of [38], which states that for all / in Hq{Q), 

wfL^m « ii^3/iii(^) + imwhn) + 11^1^3/111.(0) • 

Its proof, however, is entirely different from that of Theorem A. 2 or the inequalities 
in [36, 40] described above. (A 3D version of it can, however, be obtained using an 
argument somewhat along the lines of the proof of Theorem A. 2.) 



Appendix B. Special boundary conditions 



For the Navier-Stokes equations in 2D, when a = k, the Navier boundary conditions 
reduce to the conditions, u ■ n = uj{u) = ([10, 28, 23]).^ Here, k is the curvature of 
the boundary of a planar bounded domain and u{u) is the scalar curl of u. The natural 
extension of this observation to 3D is Lemma B.l, which involves the shape operator,'^ 



dn 

■= IT 

ov 



the directional derivative of n in the direction, f , for any vector, f , in the tangent plane. 

Lemma B.l. The boundary conditions in (1.5) reduce to those in (1.6) when A is the 
shape operator. 

Proof. Let A be the shape operator and let r be any unit tangent vector. Then since 
the shape operator is symmetric, we can write (1.5) as S{u^)n ■ r + At ■ = 0. But as 
in [2], 2S(M^)n ■ r = (curlti^ x n) ■ r — 2m^ ■ so (1.5) becomes 



(curlw"^ X n) 



dn 



At 



■ u 



0. 



□ 



In these references, the relation is written a = 2k, since 2S{u) rather than S{u) is used in the (2D 
version of) (1.3). 

''when an inward unit normal convention is used, the expression for A contains a negative sign. 
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Appendix C. Some lemmas 
We assume that is a bounded domain in having a Lipschitz boundary, F. 
Lemma C.l. Let f he a divergence- free vector field in H'^{Q) that satisfies 

S(/)n = $ on r 

in the sense of a trace, where $ lies in H^^'^{T). Then, for any vector field, g, in H^, 
we have 

- [ Af-gdx = 2 [ Sif)-Sig)dx-2 [ ^-gdS, 
Jn Jn Jr 

where A - B = Z]i<jj<3^ii% for matrices A = {aij)i<ij<3 and B = (%)i<jj<3- 
We recall the following classical lemma: 

Lemma C.2. Let u be a divergence-free vector field, of class H^{Q)^, in a bounded 
domain, Q C M.^, with a C^-boundary, T. Then, if the normal component of f vanishes 
on V , we have 

1 1 

for a constant kq depending on the domain. 
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